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Abstract. We present the first INTEGRAL results on Cyg X-3 from the PV phase observations of the Cygnus region. The source 
was clearly detected by the JEM-X, ISGRI and SPI detectors. The INTEGRAL observations were supported by simultaneous 
pointed RXTE observations. Their lightcurves folded over the 4.8 hour binary period are compatible with the mean RXTE/ ASM 
and CGRO/BATSE light curves. We fit our broad-band X-ray/y-ray spectra with a physical model, which represents the first 
such published model for Cyg X-3. The main physical processes in the source are thermal Comptonization and Compton 
reflection with parameters similar to those found for black-hole binaries at high Eddington rates. 
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1. Introduction 

The bright X-ray binary Cyg X-3 was discovered during an 
early rocket flight already in 1966 (Giacconi et al. 1967) but it 
remains still poorly understood. It is a massive system with the 
donor star and the compact object orbiting each other in a tight 
orbit. The system is embedded in a dense wind from the donor 
star, presumably a massive nitrogen-rich Wolf-Rayet star with 
huge mass loss (van Keerkwijk et al. 1992). The nature of the 
compact object is not known but recent mass estimates suggest 
it might be a black hole (e.g. Schmutz, Geballe & Schild 1996). 

The system has been observed throughout a wide range of 
the electromagnetic spectrum (e.g. McCollough et al. 1999). 
It is one of the brightest Galactic X-ray sources, displaying 
high and low states and rapid variability in X-rays. It is also 
the strongest radio source among X-ray binaries, and shows 
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both huge radio outbursts and relativistic jets. The most strik- 
ing feature in the lightcurve is a 4.8-hr quasi-sinusoidal mod- 
ulation, present both in X-rays and infrared. The modulation 
is believed to reflect the orbital motion of the binary with the 
emission from the X-ray source being scattered by the wind 
from the companion. 



2. Observations and data analysis 

On 2002 Dec. 22-23, Cyg X-3 was observed by all the 
X/y-ray instruments aboard INTEGRAL - the JEM-X (Lund 
et al. 2003), IBIS/ISGRI (Lebrun et al. 2003), IBIS/PICsIT 
(Di Cocco et al. 2003) and SPI (Vedrenne et al. 2003). 
The INTEGRAL observations were supported by simultaneous 
RXTE/PCA and HEXTE observations making possible a com- 
parison of the results from the two X-ray telescopes. At the 
time of the INTEGRAL observations, Cyg X-3 was in a rela- 
tively high state with the X-ray flux varying between 130 and 
330 mCrab, according to the RXTE/ ASM dwell-by-dwell data. 
Cyg X-3 was also observed in radio by the RATAN and Ryle 
telescopes. The results of the radio observations will be pre- 
sented elsewhere (Hjalmarsdotter et al., in preparation). 
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Fig. 1. The ISGRI 20-40 keV (top) and 40-100 keV (bottom) 
mosaic images showing Cyg X-l (left), Cyg X-3 (center) and 
SAXJ2103.5+4545 (right). 




Fig. 2. The SPI 34-52 ke V image. The brightest source in the 
lower left corner is Cyg X- 1 . Cyg X-3 is the source to the upper 
right. 



2.1. JEM-X 

The observations were performed with the second (JEM-X2) of 
the two identical X-ray monitors. The data were taken from 43 
pointings (science windows), 2200 s each, performed between 
LTD 1085.6-1087.5 (UD=fractional number of days since 2000 
Jan. 1 00.00 UT, which corresponds to MJD = LTD + 51544) 
during the 23rd rev. of INTEGRAL (2002 Dec. 22-23), which 
corresponds to 10 binary periods and the net exposure of 
~90 ks. In a half of the pointings, Cyg X-3 was in the fully 
coded field-of-view (FOV, within 2.4° from the FOV center), 
while in the rest, it was in the partially coded FOV (within 5° 
from the FOV center). In all the pointings, Cyg X-l and SAX 
J2 103. 5 +4545 (an X-ray pulsar, discovered by Hulleman, in 
't Zand & Heise 1998) were > 10° from the FOV center and 
hence the data were not contaminated by those sources. The 
offsets were rather uniformly distributed as a function of the 
binary phase. 

Source spectra were extracted individually per pointing. 
Then the average spectrum was obtained from the sum of the 
individual spectra weighted by the exposure time. The spec- 
tral response was Crab-corrected appropriately for this time pe- 
riod (instance 0004). The spectra were implicitly background- 
subtracted by a deconvolution algorithm assuming a spatially 
flat background. We used the energy range of 2.6-27 keV for 
spectral fitting. 

2.2. IBIS/ISGRI 

The ISGRI fully coded FOV is about 9 x 9°, while the partially 
coded FOV extends up to 29 x29°. Standard spectral extraction 
is at present feasible only in the fully coded FOV, therefore 
we limited ourselves to the science windows where Cyg X-3 
was at a distance of < 4.5° from the pointing direction. The 
selected ISGRI data contain about 40 science windows (of an 
average duration of 2200 s each with exception of two with 500 
s duration each) for a total duration of ~79 ks. 

Imaging analysis was performed using the current version 
of the Offline Scientific Analysis (OSA) software, using the 
procedure described in Goldwurm et al. (2003). Cyg X-3 was 
detected at a high signal-to-noise in the 15^0 and 40-100 keV 
energy bands in this as well as in previous ISGRI observations 
(Goldoni et al. 2003). The source position was obtained with 
an offset of < V with respect to the catalog position. Mosaic 
images in the 20-40 and 40-100 keV bands are shown in Fig. 
1, showing Cyg X-l, Cyg X-3, and SAX J2103.5+4545. 

Spectral extraction was performed independently for every 
science window in 24 channels linearly rebinned in the 13- 
200 keV range from a 2048-channel response matrix developed 
at CEA/Saclay. We took the source position as obtained from 
the imaging procedure and then fitted source and background 
fluxes in each energy band. The resulting individual spectra 
were added to obtain the total spectrum. A 5% systematic error 
was then added in quadrature to each channel. 

Cyg X-3 was not detected by the IBIS/PICsIT (Foschini, 
private comm.), which becomes efficient only at energies <: 250 
keV (Di Cocco et al. 2003). Given the low flux from the source 
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Fig. 3. The JEM-X 6-15 keV (red diamonds) and ISGRI 20-40 keV (blue squares) lightcurves folded over the orbital period. We 
also show the RXTEfPCA 3-15 keV data (black asterisks) from simultaneous observations as well as the RXTE/ ASM 1 .5-12 keV 
(dashed red curve) and CGRO/BAlTSE 20-100 keV (blue curve) phase dependences averaged over several years of monitoring. 
The count rates are normalized to the respective maxima. 



above 250 keV, a significantly longer exposure would be re- 
quired for detection. 

2.3. SPI 

Out of the 95 dithering pointings taken during the rev. 23 on 
the Cygnus field, 10 had to be excluded from the SPI analy- 
sis as they were either affected by strong solar activity or by 
the radiation belts. This left 85 dithering pointings with a total 
exposure of 169 ks for the present analysis. As the SPI data 
are background-dominated, a careful background substraction 
is essential in order to get reliable results, especially for weak 
sources. A time-dependent background model has been applied 
to the data, based on the saturated events seen by the detec- 
tor. The image reconstruction was done using the SPI Iterative 
Removal Of Sources program (SPIROS; Skinner & Connell 
2003). To get precise flux values, the source positions of the 
two brightest sources in the field (Cyg X-l and Cyg X-3) have 
been fixed to their catalogue values. No source confusion is 
expected in the SPI data as there are no other sources visible 
within 3° around Cyg X-3. The SPI image is shown in Fig. 2. 

For spectral extraction, 20 logarithmic bins in the 20-300 
keV energy range have been used. The instrumental response 
function used here has been derived from on-the-ground cali- 
bration and then corrected based on the Crab calibration obser- 
vation. A 5% systematic error has been added to the spectrum. 



2.4. PCA/HEXTE 

The RXTE data overlapping with the INTEGRAL observations 
are from two pointings on 2002 Dec. 22-23, of duration of 
8224 s (data set 1) and 9584 s (data sets 2 and 3), respectively. 
For the second pointing, there was a change in the number of 
PCUs used, which required breaking it into two parts, with the 
lightcurves corrected to the five PCUs. The data sets 1, 2 and 3 
then cover the binary phases of (0.86-0.06, 0.16-0.34), 0.47- 
0.67 and 0.82-1.0, respectively. Hence, data sets 2 and 3 are 
from around the maximum and minimum phases, respectively. 
The INTEGRAL spectra, accumulated from all phases, should 
have flux levels in the middle of that from these two sets. 

A 1 % systematic error has been added to the PC A spectra. 
The relative normalization of each of the data sets from the 
two HEXTE clusters with respect to the PCA data has been 
allowed free in the fits. The response matrices and background 
files have been obtained using the FTOOLS v. 5.2. 

3. Lightcurves 

From the JEM-X observations, lightcurves in four energy 
bands (3-6, 6-10, 10-15 and 15-35 keV) were created. The 
ISGRI lightcurve was extracted in the 20^1-0 keV band. The 
lightcurves were folded using the latest published ephemeris 
for Cyg X-3 (Singh et al. 2002). 

The results are plotted in Fig. 3. Two middle JEM-X bands 
were used since calibration for the 3-6 and 15-35 keV bands 
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has not been yet consolidated. For comparison, we also plot 
the folded RXTE/ ASM (2-12 keV) and CGRO/BATSE (20-100 
keV) lightcurves from monitoring of the source during 1996- 
2002 and 1991-2000 respectively. The JEM-X 6-15 keV band 
shows good agreement with the ASM data and the ISGRI 20- 
40 keV follows the shape of the BATSE curve. Modelling of 
the lightcurves will be presented elsewhere (Hjalmarsdotter et 
al., in preparation). 

4. Broad-band spectral modelling 

The /ttT£/PCA-HEXTE, JEM-X, ISGRI and SPI data were 
fitted using the XSPEC package (Arnaud 1996). We analyze 
here three RXTE spectra and the average INTEGRAL spectrum. 

We interpret the intrinsic spectra of Cyg X-3 in terms of 
Comptonization of soft X-ray seed photons, assumed here to be 
a blackbody with a temperature, T s . We use a Comptonization 
model by Coppi (1992, 1999), eqpair, described in detail 
by Gierlihski et al. (1999). This model was also used to fit 
X-ray spectra of GRS 1915+105 and Cyg X-l by Vilhu et 
al. (1999), Zdziarski et al. (2001) and by Poutanen & Coppi 
(1998), Gierlihski et al. (1999), Zdziarski et al. (2002b), respec- 
tively. In general, the electron distribution in this model can 
be purely thermal or hybrid, i.e., Maxwellian at low energies 
and non-thermal at high energies, if an acceleration process is 
present. This distribution, including the electron temperature, 
T, is calculated self-consistently from the assumed form of the 
acceleration (if present) and from the luminosities correspond- 
ing to the plasma heating rate, L^, and to the seed photons ir- 
radiating the cloud, L s . The plasma optical depth, t, includes 
a contribution from e ± pairs. The importance of pairs depends 
on the ratio of the luminosity to the characteristic size, r, which 
is usually expressed in dimensionless form as the compactness 
parameter, I = Lcrj/(rm e c 3 ), where cr x is the Thomson cross 
section and m e is the electron mass. Hereafter, the indices of ( 
have the same meaning as those of L. 

We find all the studied spectra compatible with the hot 
plasma being completely thermal, and with kT <sc 511 
keV. Then, the e ± pair production is negligible, and the ab- 
solute value of the compactness is only weakly important. 
Accordingly, we assume a constant { s = 10 (which is typical 
for accreting X-ray sources, e.g., Gierlihski et al. 1999). 

A complex issue in Cyg X-3 is the structure of its X-ray 
absorber and the presence of additional spectral components 
in soft X-rays (e.g. Molnar & Mauche 1986; Nakamura et al. 
1993). Given that our data cover energies ^ 3 keV only, we 
neglect any additional soft X-ray components and we use a 
relatively simple model of the absorber. Namely, we assume 
that an absorbing medium with the column density, Nh,o, fully 
covers the source, and another medium with the column, Nu,\, 
covers a fraction, f\, of the source. A similar model is often 
used to model the similarly complex absorption of the Seyfert 
galaxy NGC 4151 (see, e.g. Zdziarski et al. 2002a). We as- 
sume the elemental abundances of Anders & Ebihara (1982). 
We include Compton reflection (Magdziarz & Zdziarski 1995), 
parametrized by an effective solid angle subtended by the re- 
flector as seen from the hot plasma, Q, and assuming an in- 
clination of 60°. We also include an Fe Ka fluorescent line, 



which we model as a Gaussian (with the physical and equiva- 
lent widths of cr Ka and Wkc, respectively, and the peak energy 
at E^a). We allow the reflecting medium to be ionized, using 
the ionization calculations of Done et al. (1992). We define the 
ionizing parameter as £ = 4-nF\ on /n, where F[ on is the ionizing 
flux and n is the reflector density. Given the simplified treat- 
ment of the ionized reflection of Done et al. (1992), we impose 
a condition of £ < 10 4 . We assume the temperature of the re- 
flecting medium of 10 6 K. 

Given the above approximated treatment of the absorption 
and ionized reflection, the full description of the part of the 
spectrum ;S 10 keV is likely to be more complex than that given 
by our model. However, it provides a statistically satisfactory 
description of the data, including the absorbed part of the spec- 
trum, and allows us to calculate the broad-band intrinsic spec- 
trum of the source. 

During our fits, we have found that the data, covering only 
photon energies ^ 3 keV, rather poorly constrain the temper- 
ature of the seed blackbody photons. For example, we get 0.3 
keV < lcT s < 0.5 keV within 90% confidence for the RXTE data 
set 1. For simplicity, we fix it at the respective best-fit value 
for each spectrum when determining the confidence regions of 
other parameters. 

On the other hand, our data yield accurate spectral infor- 
mation only for energies ^100 keV. Thus, they poorly con- 
strain possible electron acceleration, which can be present in 
the plasma in addition to the thermal heating. Nonthermal pro- 
cesses are, in fact, clearly observed in some other spectral states 
of Cyg X-3 observed by RXTE and the CGRO/OSSE (work in 
preparation). However, allowing for the presence of nonther- 
mal electrons improves the fit only weakly, e.g. by A^- 2 ^ -2 
for the RXTE data set 1, and thus it is not required in our mod- 
els statistically. The fraction of the total power supplied to the 
plasma in electron acceleration is constrained to ^ 0.5, and the 
power law index, F acc , of the acceleration process is not con- 
strained at all at 90% confidence (typical obtained values are 
T acc ~ 2-4). We also note that the presence of nonthermal pro- 
cesses in a very similar state of GRS 1915+105 is required only 
by the data at £ 100 keV (see Fig. 3a in Zdziarski et al. 2001). 

The Fe Ka line is found to be narrow in all the RXTE spec- 
tra, with the width much below the instrumental resolution of 
the PC A. The plasma parameters obtained are given in Table 1, 
and the spectra for the 3 data sets are shown in Fig. 4. 

The INTEGRAL spectrum has been found to be rather sim- 
ilar in shape to the RXTE ones. However, unlike the PCA, the 
JEM-X data appear to require the Fe line to be broadened, with 
the corresponding decrease of A^- 2 = -12. Thus, we have de- 
cided to allow for the broadening, but, given the limited reso- 
lution of JEM-X, we kept it then frozen at the best-fit value of 
o"Kd — 0.25 keV. The resulting parameters are given in Table 1, 
and the spectrum is shown in Fig. 4. We see that the although 
the fit parameters are similar to those of the RXTE fits, the 
normalization of the JEM-X spectrum is lower than that the 
avarage of the PCA spectra by a factor of ~ 2, which is due to 
instrumental effects, see below. 
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Table 1. Model parameters 3 for the RXTE and INTEGRAL spectra. The RXTE data sets 1, 2 and 3 are from binary phases 
(0.86-0.06, 0.16-0.34), 0.47-0.67 and 0.82-1.0, respectively. The INTEGRAL spectra were accumulated from all phases. 



Data Mio N HJ fi kT s 4/A t kT h Sl/2n f E Ka W Ka E bol A X 2 h 



10 22 cm 2 10 22 cm 2 keV keV erg cms 1 keV eV erg cm 2 s 1 



RXTE(l) 


11.8^ 


256% 


61 +0 02 

u - ul -0.02 


0.37f 


71+0.02 


23 +0 03 

"•■"-0.02 


69 


1 o +0 - 2 


10000-sooo 


6 56 +0 - 04 


400! 3 " 


8.2 x 10" 9 


239/250 


RXTE(2) 




302% 


63 +0 02 


0.41f 


14 +002 

u,1 ^-0,01 


16 +002 

iu -0.02 


71 


1 o +0 - 2 


10000-5000 


6 56 +0 - 06 

U - JU _0.02 


310 +3 ° 


13.5 x 10-' 


250/249 


RXTE(3) 




330% 


63 +0 01 
"■ u -'-o.oi 


0.45f 


17 +002 


19 +003 

"■ i7 -0.02 


69 


1 4 +03 


QQQQ+1000 


6.58 + 5J;o 3 




7.0 x 10" 9 


174/226 


INTEGRAL 


16.2™ 


334+36 


44+0.02 

"•^ -0.02 


0.38f 


n 1 s+o- 05 


19 +003 

"■ i7 -0.02 


75 


1 l +0 - 2 


10000_ 3000 


fi co+0.04 


230% 


4.3 x 10- 9 


291/186 



a The uncertainties are for 90% confidence, i.e., A^ 2 = 2.71; 'f denotes a parameter fixed in the fit. 

b Calculated from the energy balance, i.e., not a free fit parameter. 

c Assumed < 10 4 in the fits. 

d The bolometric flux of the absorbed model spectrum. 




E [keV] 

Fig. 4. Deconvoled spectra of Cyg X-3. The red, magenta and cyan spectra correspond to the PCA/HEXTE data set 1, 2, 3, 
respectively. The HEXTE spectra have been renormalized to the flux level implied by the PCA. The model spectrum is shown 
only for the middle spectrum (data set 1) for clarity. The green, blue and black spectra are from the JEM-X, ISGRI and SPI, 
respectively, with the model spectrum shown in the black curve. The ISGRI and SPI spectra have been renormalized to the JEM- 
X data. The intensity levels of the INTEGRAL spectra, which were accumulated over several binary phases, should be close to 
the mean of the RXTE data sets 2 and 3. The proportions of this figure correspond to equal length per decade on each axis. 
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Table 2. The values of the relative normalizations between the 
flux level from different detectors implied by our data. 



Instruments 1 


I1-E2 


[keV] a 


ratio b 


HEXTE0/PCA(2) 


18- 


•120 


0.71 


±0.04 


HEXTE1/PCA(2) 


18- 


-120 


0.70 


±0.05 


JEM-X/PCA(2) C 


2.6 


-27 


0.33 


±0.01 


JEM-X/PCA(3) C 


2.6 


-27 


0.65 


±0.01 


ISGRI/JEM-X 


35- 


•220 


4.7 


±0.2 


SPI/JEM-X 


23- 


•100 


2 3 


±0.1 



a The energy range of the first of the two compared instruments 

used to derive the ratio. 
b The uncertainties are l-cr. 

c The JEM-X spectrum should have about the same level as an av- 
erage of the PCA (2) and (3) spectra, which implies the renormal- 
ization factor of -0.44. 



5. Comparison between individual INTEGRAL and 
RXTE spectra 

The INTEGRAL spectra were accumulated over the entire bi- 
nary orbit, unlike the three RXTE data sets (see Fig. 3). The 
mean flux levels of the INTEGRAL spectra should then corre- 
spond closely to the mean of the (extreme) RXTE data sets 2 
and 3. However, as we see in Fig. 4, the flux level of the JEM- 
X spectrum is about a half of that. Furthermore, the normal- 
izations of the HEXTE, ISGRI and SPI differ from each other, 
with the ISGRI spectrum having the flux level about wice of 
that of the average PCA. Those differences are instrumental. 
To facilitate appropriate corrections to the fluxes from various 
instruments, we list the relative normalizations between the in- 
struments in Table 2. The relative ratio between the HEXTE 
and PCA spectra is consistent with previous results. The coeffi- 
cients involving the INTEGRAL data correspond to the calibra- 
tion as of 2003 June. 

The slopes of various spectra agree with each other well, as 
shown in Fig. 4. The ISGRI and SPI data fitted by a power law 
over the energy ranges given in Table 2 yield the photon index 
of T = 3.6 + 0.1 with^ 2 /v = 34/22, 9/9, respectively. On the 
other hand, those data show hardenings at lower energies (not 
shown in Fig. 4), which appear related to residual inaccuracies 
of the present response matrices. 

6. Conclusions 

Cyg X-3 was clearly detected by all three X/y-ray instru- 
ments on board INTEGRAL. The JEM-X 6-15 keV lightcurve 
folded over the orbital period shows good agreement with 
the RXTEfPCA 3-15 keV and RXTE/ ASM 1.5-12 keV 
lightcurves. The ISGRI 20^10 keV folded lightcurve matches 
the CGR0/BA7SE 20-100 keV one. A difference in the light 
curve profile between energies above and below ~15 keV is 
indicated. 

For the first time, we fit Cyg X-3 X/y-ray spectra with a 
physical model. The main radiative processes implied by the 
INTEGRAL and RXTE data are thermal Comptonization and 



Compton reflection. The obtained intrinsic spectrum appears 
similar to that of GRS 1915+105 at a similar Eddington ratio. 

At the time of writing, there are apparent differences in 
the normalizations between the RXTE/PCA, JEM-X, SPI and 
ISGRI spectra. Calibrations and responses at this stage are con- 
stantly being improved. Cyg X-3 will be further observed with 
INTEGRAL as a part of the Core as well as the Guest Observer 
programmes. 
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